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Spontaneous breaking of time-reversal symmetry in superconductors with the -h ipy symmetry 
of the order parameter allows for a class of effects which are analogous to the anomalous Hall effect 
in ferromagnets. These effects exist below the critical temperature, T < T^. We develop a kinetic 
theory of such effects. In particular, we consider anomalous Hall thermal conductivity, the polar 
Kerr effect, the anomalous Hall effect, and the anomalous photo- and acousto-galvanic effects. 


Introduction: One of the leading candidates for p-wave 
pairing in electronic systems is Sr2Ru04. There are nu¬ 
merous pieces of experimental evidence that the super¬ 
conducting state of Sr2Ru04 has odd parity, breaks time 
reversal symmetry and is spin triplet [1-6] . An order pa¬ 
rameter consistent with these experiments is given by the 
chiral p-wave state [7] which is an analog of ^He-A. It has 
the form Aq,^ (p) ~ ± ipy where Aq,/? (p) is the Fourier 

transform of AQ,/3(r — r'). However, the observation of 
power laws in specihc heat [8] and NMR [9] , the absence 
of electric currents along edges [10], and the absence of 
a split transition in the presence of an in-plane magnetic 
field [11] are inconsistent with the theoretically expected 
properties of a simple chiral superconductor. Consider¬ 
ation of additional experimental manifestations of spon¬ 
taneous breaking of time-reversal symmetry in px + ipy 
superconductors may clarify the nature of superconduct¬ 
ing state in Sr2Ru04. 

Due to spontaneous breaking of time-reversal sym¬ 
metry, px + ipy superconductors must exhibit anoma¬ 
lous transport phenomena similar to those which exist 
in metallic ferromagnets (see Refs. [12] and [13] for a re¬ 
view) . In this article we develop a theory of several such 
effects in px + ipy superconductors: the anomalous Hall 
effect, polar Kerr effect for microwave radiation, anoma¬ 
lous Hall thermal conductivity, and anomalous photo- 
and acousto-galvanic effects. 

It should be noted that p-wave superconductivity ex¬ 
ists only in the clean regime, I > ^, where electron 
transport may be described semiclassically. Generally, in 
the semiclassical regime there are three contributions to 
anomalous transport phenomena: skew scattering, side 
jumps, and the intrinsic contribution. The side jump 
contribution arises from the shift of the center of mass of 
electron wave packets during the scattering events, while 
the intrinsic contribution is related to the anomalous ve¬ 
locity due to Berry curvature. The magnitude of these 
contributions is independent of the mean free path. In 
contrast, the magnitude of the skew scattering contribu¬ 
tion is proportional to the quasiparticle mean free path 
1. As a result, the skew scattering contribution exceeds 
the intrinsic and side jump contributions by a factor 

PfI 1- Here p_F is the Fermi momentum. There¬ 
fore in this article we will take into consideration only 


the skew scattering contribution. We focus on anoma¬ 
lous transport phenomena in the vicinity of the critical 
temperature, where quasiparticles play a major role. 

Kinetic scheme: Transport theory in conventional time 
reversal invariant superconductors was developed long 
ago (see for example reviews Refs. [14] and [15]). Below 
we generalize this approach to superconductors without 
time reversal symmetry, which exhibit anomalous trans¬ 
port phenomena. In the clean regime, I ^ and at suf¬ 
ficiently low frequencies, a; <C |A|, where |A| is the mod¬ 
ulus of the order parameter, the quasiparticle dynamics 
can be described by the Boltzmann kinetic equation for 
quasiparticles 


dnp{r,t) dip dup dip dnp 


dt 


dp dr dr dp 


= L 


( 1 ) 


where 


+ v-p^, ep = i7C2-f|A| 


(2a) 


fp - Cp + $ + , ?p - ep. (2b) 

In Eq. (2) m is the electron mass, while Ps and $, are 
given by 

= (I> =^dtx + e(l), (3) 

2 c 2 


where x is the order parameter phase, and (j) and A are 
the scalar and vector potentials. From Eq. (3) one ob¬ 
tains the equation for the acceleration of the condensate. 


5tPs = eE-k V$. (4) 

Equations (l)-(3) should be supplemented by the expres¬ 
sion for the current density, 

eN f 

j = -Ps + e / dFvrip, (5) 

m J ^ 

and by the charge neutrality condition, 

= J dr ^ Up, ( 6 ) 

that relates the gauge invariant scalar potential and the 
odd in ^ part of the quasiparticle distribution function. 
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the self-consistency equation for the order parameter. 
Here dT = Vcfip/{2TTfi)^ {V is the volume of the sam¬ 
ple) and V = d^p/dp. 

We work in linear response to external perturbations, 
and neglect corrections to equilibrium value of |A|. The 
collision integral 1st = + lit Eq. (1) describes 

both elastic and inelastic scattering. We will assume that 
Te ^ r, where and r are inelastic and elastic mean free 
time respectively. Therefore the main contribution to 
the aforementioned anomalous effects comes from elastic 
scattering, which is described by the collision integral 

1st = J {Wpp,np, - S {ip - ip,) dT'. (7) 

Skew scattering of quasiparticles corresponds to the 
part of scattering probability in Eq. (7) that is associ¬ 
ated with breaking of time reversal symmetry, SWpp, = 
Wpp! — W-pi-p 7 ^ 0. Thus, all the aforementioned ef¬ 
fects are proportional to SWpp,. Skew scattering arises 
beyond the lowest Born approximation for the scattering 
amplitude. Below we consider point-like impurities. In 
the normal state such impurities scatter electrons only 
in the s-wave channel and do not cause skew scattering. 
Therefore in the superconducting state skew scattering 
of quasiparticles is entirely due to the breaking of time 
reversal symmetry by the Px + Wy order parameter. The 
elastic scattering probability for quasiparticles with en¬ 
ergy e can be characterized by ^ = ^p, ^' = ^p, = and 
the asimuthal angles p, p', which define the direction 
of p and p' in the cry-plane. For simplicity, we assume 
cylindrical Fermi surface and obtain for the scattering 
probability (see appedix for details). 

VFpp' = Wo + Will- cos {ip-v^' + 26,)]. (8) 

Here 6, is the energy-dependent scattering phase shift. It 
is related to the s-wave scattering phase shift (5„ in the 
normal state by 


J^=arctan— 

- A2 


(9) 


We assume weak impurities, for which Sn ~ tan(5„ = 
—TTiyVo is small. Here is the density of states on the 
Fermi level and Vq is the impurity pseudo-potential [16]. 
In this case Wq and Wi are given by 


W^o(C,r) 


C(e) (g + g')^ 

2z/r 2e^ 


2iyT 


(10a) 

(10b) 


Here = 2TrniuVo, with rii being the impurity density, 
is the elastic scattering rate in the normal state. The co¬ 
efficient ^(e) = (e^ — A^)/[e^(l-1-5^) — A^] represents the 
enhancement factor of the quasiparticle scattering cross- 
section over the normal state value. The first term in 


Eq. (8), Wo given by Eq. (10a) has the same structure 
as in s-wave superconductors. It describes scattering 
only within the same (particle-like, ^ > 0, or hole-like, 
^ < 0) branch and does not lead to branch imbalance 
relaxation. The second term, Wi in Eq. (8) is absent in 
s-wave superconductors. It leads leads to both skew scat¬ 
tering and scattering between branches of quasiparticle 
spectrum with different signs of The skew scattering 
cross-section, described by the sin((p — p') sin 26, term in 
Eq. (8), is energy-dependent. It follows from Eqs. (8), 
(9), and (10b) that it changes sign when impurity poten¬ 
tial Vo changes from repulsive to attractive. 

Below we consider linear response to several external 
perturbations and look for the quasiparticle distribution 
function in the form Up = -I- np\ where is a 

locally equilibrium Fermi distribution, and rip'^ describes 
the deviation from equilibrium. Noting that the collision 
integral (7) is nullified by an arbitrary function n^^^ip) 
we write the linearized Boltzmann equation in the form 

5(p) = J dT'WppKnW - nW)5(ep - e^,), (11) 

where the specific form of the source 5(p) depends on 
the type of perturbation. 

Anomalous Hall thermal conductivity: We first con¬ 
sider the Hall component of the thermal conductivity K^y 
which describes the heat flux perpendicular to the direc¬ 
tion (x-axis) of the temperature gradient. In this case 
the source term in Eq. (11) has the form, 

5(p) = -|;v.vr—. (12) 

The expression for the heat flux is 

jQ = J dTep^n<^l (13) 

Note that dcp/dp = v^/e is the group velocity of the 
quasiparticles while v is the bare velocity as in a normal 
metal, jvj = vp- The solution of Eqs. (11), (12) has the 
form 

«p = [as(e)sin(y-bac(e) cos(y]. 

The Hall component of the thermal conductivity ten¬ 
sor, Kxy, is determined by as(e) in the above expression, 
which is given by as{e) = a(e)/[6^(e) -I- a^(e)], with 

a(e) = ^ sin2h„ (14a) 

For weak impurities, |(5„| <C 1, we obtain, close to T, 

/ A V 

f^xy — ^ J 6n: (Ih) 
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where k = ir'^i/TD/S (with D = v\tI2 being the diffu¬ 
sion constant) is the normal state thermal conductivity. 

Polar Kerr effect: Next we consider a linearly polar¬ 
ized electromagnetic wave at normal incidence to the xy 
surface of Px + ipy superconductor. The reflected wave 
is elliptically polarized with the major axis rotated with 
respect to the incident one by the polar Kerr angle [17] 

{1 — n'^ + k‘^)Ak + 2nKAn 

^ ^ (1 - n2-h k2)2 + (2„k)2 ’ ^ > 


where n and k are, respectively, the real and imaginary 
part of the refraction index and 


An -|- iAk = 


Itt (n — iK)axy 
LO rp + ’ 


(17) 


where CFxy is the complex ac conductivity. 

In this case the electric field is uniform in the direction 
parallel to the surface of the sample, d) = 0, and the value 
of Ps(t) is determined by Eq. (4). The diagonal part of 
the conductivity is given by [14] 


^XX 


ctd + 


iNs{T) 

cu 


(18) 


where Ns{T) is the temperature dependent superfluid 
density and ctd = e^uD is the Drude conductivity. In 
contrast to the thermal conductivity consideration, in the 
present case n^°^ = l/(exp[ep -I- v • Ps{t)]/T + 1) gives a 
nonvanishing contribution to the current response (via 
Ns{T)) because the superfluid momentum depends on 
the electric field. The Kerr angle 9 is determined by the 
value of the Hall component of conductivity axy 

To find axy we seek the solution of Eq. (1) in the form 
n = n (o)(ep/T) -I- rip ^ The source in Eq. (11) becomes 

5(p) =- ev • E-^, (19) 

where the external electric field E is along x-direction. 
The nonequilibrium distribution rip^^ has the form 


ffi) = - 


= -eEvp- 




[,5s(e)sin(/j-|-/3c(e)cos(^]. (20) 


The Hall conductivity depends only on the function (e), 
which is given by ffi^e) = a(e)/{[6(e) —ja;]^-|-a^(e)}, with 
a(e) and 6(e) being defined in Eq. (14). Substituting 
Eq. (20) into Eq. (5) we obtain the Hall conductivity in 
the weak impurity limit, [ttz/VqI ^ 1, in the form 

■»( ) “ cosh“(VJ^A/2r) 

X-- . ( 21 ) 

{—iujTxffx'^ + 1 -I- x2 -I- 3/4)2 

where x = [■Cl/^- At temperature close to Tc and at low 
frequencies, lot <C 1, this expression yields 


This result was derived assuming px + ipy symmetry of 
the order parameter. In the Px — ipy state the Hall con¬ 
ductivity axy has opposite sign. It also changes sign if 
the impurity potential Vq changes from repulsive, < 0, 
to attractive, > 0, in agreement with Ref. [18]. Note 
that our result for the low frequency Hall conductivity, 
Eq. (22), is proportional to the elastic mean free time t 
and to the density of quasiparticles. 

There is another contribution to axy associated with 
the existence of the transverse component of the super¬ 
fluid velocity Vy ^ Px, which is proportional to the con¬ 
densate acceleration in the x-direction. It may not be 
obtained within the present formalism that is based on 
the Boltzmann kinetic equation for the quasiparticles. At 
T ~ Tj, this contribution is smaller than the quasiparti¬ 
cle contribution, Eq. (22). However at T -^Tc when the 
quasiparticle contribution becomes exponentially small 
in Eq. (21) it becomes the dominant contribution. The 
requirement for this contribution to exist is violation of 
Galilean invariance in the system. Thus it should exist 
in any crystalline superconductors with px + ipy symme¬ 
try [19, 20]. It can also be caused by electron-impurity 
scattering. In this case this contribution is inversely pro¬ 
portional to the electron mean free time [18]. 

Hall effect for normal current injection: Let us now 
consider a normal metal/pj, -I- zp^-superconductor junc¬ 
tion, through which a steady current is flowing. At 
T A this situation was considered in Ref. [21]. In 
this regime conversion of normal current to supercurrent 
is mediated by multiple Andreev reflections. Here we 
work near the critical temperature and consider a setup, 
in which the normal current is injected into the super¬ 
conductor in the x- direction, as shown in the inset in 
Fig. 1. In this case the conversion of quasiparticle current 
to the supercurrent occurs in the superconductor. Just 
as in the case of s-wave superconductor, near Tj,, the 
electric field penetrates into superconductor to a large 
distance Lq ^ I, which is determined by the relaxation 
of imbalance between the populations of quasiparticles 
in electron-like, ^ > 0, and hole-like, ^ < 0 branches 
of the spectrum (see for example Ref. [14] and refer¬ 
ences therein). The new feature of normal current in¬ 
jection that appears in px + ipy superconductors is that 
skew scattering of quasiparticles generates nonequilib¬ 
rium current that is perpendicular to the electric field. 
Another aspect is that, in contrast to s-wave supercon¬ 
ductors, impurity scattering leads to branch imbalance 
relaxation even if the magnitude of the order parameter 
jAj is isotropic in the Fermi surface. Below we assume 
that the inelastic scattering rate is smaller than 1/r and 
thus impurity scattering gives the dominant contribution 
to branch imbalance relaxation. 


^xy 


Til A 


( 22 ) 


In linear response we write the quasiparticle dis¬ 
tribution function in the superconductor in the form 
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(cp/T)+np \ This yields the source term in Eq. (11) 


‘5(p) 



e 


clr 


(23) 


At length scales in excess of the mean free path we may 
employ the diffusive approximation. With the aid of 
Eq. (5) the Hall current jy can be expressed in the form 


jy{x) = -iei^DSn 



(24) 


where na(C, x) is the antisymmetric in ^ part of the distri¬ 
bution function averaged over the momentum directions. 
The latter satisfies the diffusion equation with relaxation 




1 


na{^,x), 


(25) 


with energy dependent relaxation rate Tq^(0 = 

-I- 2A^)/^^. The solution of Eq. (25) is 


na{^, x) = 0) exp [-x/Lq{^)] , (26) 


where Lq (^) = \/Dtq{^) is the energy-dependent branch 
imbalance relaxation length. 

We work in the vicinity of the critical temperature Tc, 
where for typical thermal quasiparticles ^ ~ T, (|^| ^ 
A), the relaxation lengths are long, Tq(^) = ^ I- 

These quasiparticles diffuse into the bulk of the supercon¬ 
ductor and contribute to the gauge-invariant potential 4) 
given by Eq. (6). The boundary value of the nonequi¬ 
librium quasiparticle population, fia{^,0) in Eq. (26) is 
obtained by matching the solution of Eq. (25) with the 
solution of diffusion equation with energy relaxation for 
the electrons in the normal metal. The result depends on 
both the inelastic mean free path in the normal metal, 
le, and the branch imbalance relaxation length Lq in the 
superconductor. For 1^ ^ Lq the boundary condition is 


»^a(?,0) = sign(^) 


eE,{0)LQ{O 

4Tcosh^(^/2T)’ 


(27) 


where (0) is the electric field in the normal metal gen¬ 
erating the steady current. Here we used the fact that 
in the stationary case E = — V$/e, which follows from 
Eq. (4). Using this relation and substituting Eqs. (27), 
(26) into Eqs. (6) and (24) we obtain the spatial distri¬ 
butions of the electric field Ex{x) and the Hall current 
jy{x) in the superconductor, 


Ex{x) 


E^{0)Fo 



(28) 


Jy{x) = aDE,{0)6r.(^^'^ (29) 

where (Lq) = 2\n2{Tl/ A) and the functions are 
defined as 


Fn{x) 


r, 

Jo ^cosh^(?/) 


exp 



(30) 


F 



FIG. 1. Plot of the functions F- 2 {x) (solid line) and Fo{x) 
(dashed line) in Eq. (30). The inset shows a schematic setup 
of the normal current injection experiment. Electric current 
is injected into the superconductor S from the normal metal 
N along along the a;-axis. Skew scattering of quasiparticles 
generates an anomalous Hall current in the y direction. 


and are plotted in Fig. 1. The spatial distributions of 
the Hall current jy{x) and the electric field Ej.{x) are 
drastically different, and cannot be related by a local 
Hall conductivity cjxy. At relatively short distances, 
I X {Lq), we see from Eq. (29) that jy{x) oc 
aDEx{0)Sn{A/T)'^{LQ)/x, so that the Hall current is 

Iy = J dxjy{x) aDE,,{0)Sn (LQ)ln^^^. 

Anomalous photo- and acousto-galvanic effects: When 
an electromagnetic or an acoustic wave propagates 
through a conductor it generates an anisotropic in mo¬ 
mentum p distribution function. The induced current 
density is proportional to the rate of the the momentum 
transfer from the wave to the electron system [15], 

1 — Try 

Here / is the rate of momentum density transfer due to 
the wave adsorption, and x is the direction of the wave 
propagation. In -h ipy superconductors an anomalous 
current in the y direction is generated. Considerations 
similar to those leading to Eq. (22) near yield 

^xy ^xx rji 

Finally, we note that all anomalous transport phe¬ 
nomena discussed above are driven by the underlying 
symmetry of the superconducting state. Therefore they 
should exist in any superconductor whose order parame¬ 
ter breaks time reversal symmetry, see for example [22- 
26] . Although our consideration focused on -|- ipy ma¬ 
terials we believe our approach is applicable to other su¬ 
perconductors with broken time-reversal invariance. 
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Appendix: Derivation of the scattering probability 

In this appendix we obtain the scattering probability 
Wpp' in Eq. (8) that gives rise to anomalous transport 
properties. For simplicity we consider point-like impuri¬ 
ties whose scattering matrix elements are independent of 
momentum, Vpp' — Vq. Skew scattering appears beyond 
the lowest Born approximation as a result of a particular 
structure of the order parameter in Px + ipy superconduc¬ 
tors. We note that even weak impurities, \vVq\ <C 1, lead 
to the existence of bound states in px -f ipy superconduc¬ 
tors with binding energies ~ IA| Ii/VbP [27]. Therefore 
low energy quasiparticles undergo resonant scattering. 
We show below that this also results in energy-dependent 
skew scattering. 

The scattering properties of a single impurity are de¬ 
scribed by the T-matrix. In a superconductor, for each 
initial and final momenta p' and p the T-matrix acquires 
an additional 2x2 structure in the Nambu-Gorkov space. 
We denote these 2x2 matrices by Tpp'. The quasiparticle 
scattering amplitude is given by the on-shell matrix ele¬ 
ment Tpp/ = (p|T|p'), with Ip) being the two component 
quasiparticle wave function. Using the Fermi Golden rule 
the scattering probability can be expressed as 

Wpp/ = 27rni|Tpp/p (A.l) 


and the Green function G'o(p) = (^E — H(p)^ 
1 


-1 


Go)?) = 


E^-el 


E + ^p 

Ae-^v’P E-^p 


(A.5) 


Replacing with V J d^p/(27rh)^ ^ f I d(pp/2TT, 

where V is the volume of the sample and v is the den¬ 
sity of state at Fermi surface, and using Eq. (A. 3) and 
Eq. (A.5), we obtain for the T-matrix, 


f(T) 


_ Yo _ 0 

l+T^vVoE/y/A'^-E'^ 

0 __ 

I-ttvVoE/^/A^-E'^ 


The pole of the T-matrix describes the subgap 
bound state [27] with the binding energy Eb = 
A/-\/l -I- For the energy in the continuum, 

iTj > A, we obtain the T-matrix in the form 


f{E) 


f{E)E^^ 

0 


0 ^ 

-/(T)e-*^- ) 


(A.6) 


where (5^; is given by Eq. (9) and 


f{E) 


^0 

V^l + (^i/Uo)2T2/(£;2 - A2) ■ 


(A.7) 


The quasiparticle scattering amplitude is given by the 
on-shell matrix elements between two quasiparticle states 


with energy e = J IS? + Q, = J IS? + , 


Tpp/ = (p|f(e)|p'). 


(A.8) 


In a Pa; + iPy superconductor the Bogoliubov amplitudes 
carry a momentum-dependent phase tpp, 


where rii is the impurity concentration. The T-matrix 
obeys the Lipmann-Schwinger equation 

Tpp/ = Upp, -f ^ Vpp^-Goip")fp„p,, (A.2) 

p" 

For point-like impurities, Upp, = Vqt^ (with fi being 
the Pauli matrices in the Nambu-Gorkov space), the T- 
matrix depends only on the energy E. Thus Eq. (A.2) 
simplifies to 

f{E) = Voh + Uofg ^ Go(p)T(T) (A.3) 

p 

The order parameter of Px + ipy superconductor can 
be expressed as A(pa; + iPy)/PF = Aexp(jipp) where ipp 
is the momentum-dependent phase such that cos tpp = 
Px/pe and sini^p = py/pp- This leads to the BCS Hamil¬ 
tonian 


H{p) 


^p \ 

-Cp ) ’ 


(A.4) 


(p| = (upe *'^'’,Up), (A.9) 

where Up and Vp have the standard BGS form and are 
independent of the momentum direction, 


Up — ^(1 + ^p/cp) /2, Up — Cp/fp) /2 (A.10) 


Substituting Eqs. (A.9) and (A. 10) into (A.8) we get 


T, 


pp' 


— f (upUp/ UpUp/) -f 2zipUp/UpUp/ 

X [1 - cos(2dE -t tpp - (/9p/)] I (A.ll) 
l-cos(25.+ i;/p-(,/p-) (A.12) 


Substituting this result into Eq. (A.l), and recalling that 
for weak impurities, [ttpI/qI ^ 1; the scattering rate in the 
normal state is = 27rniPUg^ we obtain the scattering 
probability in the form of Eq. (8) and Eq. (10). 
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